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ABSTRACT 

Persistent activity of magnetars is associated with electric discharge that continuahy injects rel- 
ativistic particles into the magnetosphere. Large active magnetic loops around magnetars must be 
filled with outflowing particles that interact with radiation via resonant scattering and spawn electron- 
positron pairs. The outflow energy is processed into copious until the plasma enters outer parts 
of the loop where the magnetic field is reduced below 10^"^ G. In the outer zone, photons scattered 
by the outflow do not convert to pairs and the outflow radiates its energy away. The escaping 
radiation forms a distinct hard X-ray peak in the magnetar spectrum. It has the following features: 
(1) Its luminosity L — 10'^^ — 10'^^ erg s~^ can easily exceed the thermal luminosity from the magnetar 
surface. (2) Its spectrum extends from 10 keV to the MeV band with a hard spectral slope, which 
depends on the object inclination to the line of sight. (3) The anisotropic hard X-ray emission exhibits 
strong pulsations as the magnetar spins. (4) The emission spectrum typically peaks around 1 MeV, 
but the peak position significantly oscillates with the spin period. (5) The emission is dominated by 
the extraordinary polarization mode at photon energies below ^ 1 MeV. (6) The decelerated pairs 
accumulate and annihilate at the top of the magnetic loop, and emit the 511-keV line with luminosity 
Lann ^ O.IL. Features (l)-(3) agree with available data; (4)-(6) can be tested by future observations. 
Subject headings: plasmas — stars: magnetic fields, neutron — X-rays 



1. INTRODUCTION 

Besides spectacular outbursts, magnetars produce per- 
sistent or decaying X-ray emission with luminosity L ^ 
10^'* — 10^^ erg s~^. Two peaks are observed in their 
X-ray spectra, with comparable luminosities. The first 
peak is near 1 keV; it is associated with thermal emission 
from the neutron star surface. The second peak is above 
100 keV. Its low-energy slope (between 10 and 100 keV) 
was observed in 7 magnetarsQ (Kuiper et al. 2008; Enoto 
et al. 2010), with a typical photon index F ~ 1 — 1.5. 
The emission shows pulsations with the rotation period; 
the pulsed fraction approaches 100% at high energies. 

The present paper focuses on the mechanism of the 
hard X-ray component. The puzzle of hard X-ray emis- 
sion was previously discussed in several works (Thomp- 
son & Beloborodov 2005; Heyl & Hernquist 2005; Be- 
loborodov & Thompson 2007; Baring & Harding 2007; 
Lyubarsky & Eichler 2008). In general, nonthermal lu- 
minosity from an isolated neutron star must be fed by 
some form of energy release in its magnetosphere. In 
magnetars, this process occurs in the closed magneto- 
sphere via electric discharge (Beloborodov & Thompson 
2007). The created relativistic particles move along the 
magnetic field lines and may radiate their kinetic energy 
in two ways: 

1. Particles moving toward the neutron star will hit 
its surface and deposit their energy in a dense layer, 
which can radiate bremsstrahlung photons. 

2. Particles can radiate their energy in the magneto- 
sphere by scattering photons around the star. 

^ Non-detections of similar hard X-ray emission in other 
magnetars are inconclusive because of insufficient sensitivity 
(W. Hermsen, private communication). 



We consider in this paper magnetic loops that extend 
to radii r ^ 5i?, where R is the radius of the neutron star 
(Section 2.1). We assume that relativistic particles are 
injected in the low parts of the loops where the magnetic 



field is very strong, B ^ Bq 
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As explained below, the particles cannot immediately ra- 
diate their energy via scattering. Instead, the energy is 
adiabatically processed into a standard outflow along the 
magnetic field lines. This processing occurs via creation 
of secondary pairs whose number grows by a factor 
of ^ 10^ as the flow moves away from the star toward 
the top of the magnetic loop (Figure 1). Pair loading 
ends where B sa 10^"^ G; at this point the outflow be- 
gins to lose energy to escaping photons. We find that 
the outflowing plasma radiates all its energy in the zone 
3 X 10^^ < B < 10^^ G, in the upper parts of the loop. 
The released radiation is practically independent of the 
details of energy injection near the star because of the 
energy processing in the pair-loading zone B > 10^^ G. 

The paper is organized as follows. Section 2 describes 
the relativistic outflow, its formation and dynamics. It 
is approximated by a simple analytical model, which is 
confirmed by detailed numerical simulations in the ac- 
companying paper (Beloborodov 2012, hereafter B12). 
Radiation emitted by the outflow is calculated in Sec- 
tion 3. It is found to have an extended hard spectrum 
between 10 keV and ~ 1 MeV and form a pronounced 
hard X-ray peak in the magnetar spectrum. Section 4 
explores the anisotropy and pulsations of radiation pro- 
duced by this mechanism. 

2. RELATIVISTIC FLOW IN THE J-BUNDLE 

2.1. Source of energy 

The radiative mechanism described in this paper is in- 
dependent of the energy source as long as the source 
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Figure 1. Sketch of an activated magnetic loop. Relativistic 
particles are injected near the star where B > Bq = 4.4 X 10^^ G. 
Large multiplicity M ~ 100 (Equation [TlJ develops in the adi- 
abatic zone B > 10^^ G (shaded in blue). The outer part of the 
loop is in the radiative zone; here the scattered photons of energy 
E = husc escape and form the hard X-ray spectrum that is cal- 
culated in Section 3. The outflow decelerates (and annihilates) at 
the top of the loop, shaded in pink; here it becomes very opaque to 
the thermal kcV photons flowing from the star. Photons reflected 
from the pink region have the best chance to be upscattered by the 
relativistic outflow in the lower parts of the loop, and control its 
deceleration (Section 2.2). 

generates relativistic particles near the magnetar where 
B > 10^^ G. We here briefly describe the expected 
source. 

Like the sun, magnetars are believed to have twisted 
magnetospheres that are deformed by surface motions 
(Thompson et al. 2000; 2002). Beloborodov (2009, here- 
after B09) developed electrodynamic theory for the dis- 
sipative, twisted magnetosphere attached to a conduct- 
ing sphere (neutron star). The theory predicts that the 
magnetic twist concentrates on field lines that extend far 
from the conducting star, forming an extended bundle 
of electric currents. This "j-bundle" is heated by ohmic 
dissipation and emits radiation. The currents are nearly 
force-free, j x B = 0, and sustained by a longitudinal 
voltage $ along the magnetic field lines. Net current cir- 
culating through the magnetosphere, /, generates ohmic 
power L = /$ that feeds the observed activity. 

Similar ohmic heating occurs in ordinary pulsars, and 
it may be useful to compare them with magnetars. In 
ordinary pulsars, the magnetospheric twist is pumped at 
the light cylinder by the rotation of the star, and ohmic 
dissipation occurs on the open field lines. The electric 
current is then roughly given by / ~ cft/Rf^ (where Ric 
is the light-cylinder radius and fl is the magnetic dipole 
moment of the star), and voltage $ can exceed 10^^ V. 
The dissipated power /$ is always smaller than the spin- 
down luminosity of the star. In magnetars, surface mo- 
tions twist the closed magnetosphere. Then the electric 
current may be as large as / ~ cfi/R^, where R is the 
star radius. The characteristic voltage is $ ^ 10^ V (Be- 
loborodov & Thompson 2007). The dissipated power /$ 
is typically much larger than the spindown power. 

Ohmic dissipation tends to gradually remove electric 
currents from the closed magnetosphere as described in 
B09. This process creates a growing "cavity" with van- 
ished current density j = 0. The currents have longest 



lifetimes on magnetic field lines with large apex radii 
^max ^ R- As a result, the magnetospheric activity 
becomes confined to the bundle of extended field lines. 
This theoretical picture agrees with observations (see Be- 
loborodov 2011 for a review). 

Consider magnetic field lines that extend sufficiently 
far from the star into the region where B is smaller than 
a given value Bi. Luminosity that can be generated on 
these field lines may be estimated as follows. Approxi- 
mating magnetic field by a dipole configuration with a 
dipole moment fl, one finds that the field lines reaching 
the region B < Bi carry the magnetic flux Fi = 27r/i/i?i 
filling the hemisphere of radius Ri = [fi/BiY^^. If the 
field lines are twisted with amplitude ■)/', they must carry 
electric current (according to Stokes theorem, see B09 
and Appendix C), 
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The twist amplitude ^ is measured in radians. It has the 
meaning of relative azimuthal position of the northern 
and southern footpoints of the magnetospheric field line. 
The generated power is L = which yields 
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Strong twists tend to inflate the magnetosphere (Wolfson 
1995; Parfrey et al. 2012a,b). As long as i/i < 1, this 
effect is modest (it scales as ip"^ a,t < 1, see B09), and 
the poloidal field may be well approximated by a dipole 
configuration with dipole moment fl. 

2.2. Interaction of outflowing particles with radiation 

field 

Discharge with voltage $ ^ 10^ V injects electrons (or 
positrons) with high Lorentz factors 7 ^ e^/nieC^ ^ 103. 
The relativistic plasma created near the star expands 
along the magnetic field and forms a relativistic out- 
flow, resembling the outflow along open field lines in or- 
dinary pulsars (however, here plasma moves along closed 
field lines and is trapped in the magnetic loops around 
the neutron star). Magnetars are hot and bright; their 
dense radiation exerts a strong drag on the magneto- 
spheric particles and controls the outfiow velocity. On 
the other hand, the plasma significantly changes the ra- 
diation field around the magnetar. The interaction of the 
flowing plasma and radiation may be described as a 
self-consistent radiative transfer problem. This problem 
is solved numerically in the accompanying paper (Be- 
loborodov, in preparation, hereafter B12). The results 
can be summarized as follows. 

The outflow interacts with radiation via resonant scat- 
tering; other processes turn out unimportant. Consider 
an outflowing electron (or positron) with Lorentz fac- 
tor 7 = (1 — /32)-i/2^ and a target photon of energy Et 
propagating at an angle "i? with respect to the outflow 
direction. Resonant scattering can occur if the photon 
energy in the electron rest frame, 



E ^-f{l- Pcos'd)Et, 
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matches Landau energy Huje, 



E = hion = hi 



The relativistic outflow in the low parts of magnetic loops 
interacts with thermal keV radiation that has been re- 
flected by plasma trapped at the top of the magnetic 
loops (Figure 1). The outflow sees the reflected photons 
as main targets because they propagate at most favorable 
angles for resonant scattering (Appendix A). 

Scattering of thermal keV radiation gives a huge energy 
boost to photons, on average by a factor ^ 7^. The rel- 
ativistic outflow experiences significant energy losses by 
scattering a tiny fraction of photons around the magne- 
tar. The outflow adjusts so that its scattering rate is just 
enough for the self- consistent deceleration. In essence, 
the relativistic outflow moves fast enough to interact with 
photons of energy E ~ 10/cT (the low-density exponen- 
tial tail of the thermal spectrum), and slow enough to 
not interact with the main peak of the thermal spec- 
trum E ~ 3A:T. This condition, together with Equa- 
tion (HI, determines that the scattering plasma moves 
with Lorentz factor 
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The self-regulation of Lorentz factor is a key feature 
of the outflow that allows one to calculate its emission 
spectrum (Section 3). This feature is illustrated in Fig- 
ure 2 that shows the dynamics of a relativistic electron 
resonantly interacting with radiation. The electron was 
injected with Lorentz factor 7 ~ 1000 at i? = 10i?Q and 
the figure shows its trajectory in the B-^ plane. As the 
electron moves outward along the magnetic field line, it 
scatters photons and loses energy. Each scattering event 
is a step down in 7 at fixed B, and each free path between 
subsequent scatterings is a step down in B at constant 7. 
In this figure, the mean free path €(&) is plotted between 
scattering events; a full Monte-Carlo simulation will be 
shown in Section 3. 

The average effect of scattering may be described as 
a drag force acting on the electron, J"(i?,7). This force 
!F is sensitive to the electron position in the B-^ plane. 
As indicated in the figure, the drag is enormous above 
the electron trajectory and negligible below the trajec- 
tory. In essence, the electron surfs the steep slope of 
J^{B, 7) so that it stays outside and around the strong- 
drag region of the B-^ plane. The relevant dimensionless 
parameter measuring the drag effect is roughly given by 
(see Appendix A), 
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where "iJmax is the maximum angle of target photons 
with respect to the electron velocity; cosi?max = —0.5 
is assumed in Figure 2. The electron surfs with D ^ 
|dln_B/(ilnr| ^ 3, which corresponds to y = 6 — 7 and 
7 « 100 5(fcr/0.5 keV)"^ This condition defines a line 
in the B-'j plane that quite accurately describes the elec- 
tron trajectory, confirming Equation ([5]). 



1000 



100 



10 



- strong drag 




0.01 



Figure 2. Deceleration of a relativistic electron moving along the 
magnetic field line away from the star and interacting with ambient 
radiation field of temperature kT 0.5 keV (see text). Each step 
down in 7 is a result of one scattering. 



2.3. Two zones of the outflow deceleration 
2.3.1. Pair loading zone 

The relativistic outfiow created in the ultra-strong 
magnetic field must decelerate as it moves to weaker 
fields, according to Equation ([SJ. However, it is initially 
unable to radiate its energy away, because the scattered 
photons quickly convert to pairs that join the outflow. 
The deceleration occurs because the flow kinetic energy 
is shared by more particles. 

The mean expectation for the energy of a resonantly 
scattered photon is given by (see Equation IA2[) 



E = "fscbnieC^ q{b), 
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which states that the photon energy in the electron 
frame, E — bnieC^ (Equation |4|), is blueshifted in the 
lab frame by the typical Doppler factor ~ 7sc. The 
factor q < 1 is a correction due to electron recoil in 
scattering. In weak fields g ~ 1, and in strong fields 
(6^ 1) q <C 1, consistent with the evident re- 

quirement E < 7scTOeC^- Combining Equation ^ with 
Equation ^ , one finds that photons capable of convert- 
ing to pairs, E > 2meC^, are generated in the region 
of rather strong magnetic field b > bo ^ 1/4. A more 
accurate formula for bo is given in Appendix B. 

Photons scattered in the zone b > bo do not escape. 
Their fate depends on the polarization state: || photons 
directly convert to , and _L photons split in two daugh- 
ter II photons (see Appendix A for the definition of || and 
_L polarization states). Conversion occurs immediately 
once the || photon satisfies the threshold condition. 
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where ■(? is the photon angle with respect to the magnetic 
field. The angle is initially small after scattering, ~ 
7~^, however it quickly grows as the photon propagates 
through the curved magnetic field. Therefore, almost all 
II photons of energy E ^ 2mf.<? produced in the zone 
b > bo are quickly absorbed. This fact is observed in our 
numerical simulations and discussed in Appendix B (cf. 
Medin & Lai [2010], where a similar problem was studied 
for the polar cap near the magnetic axis). 

In summary, the outfiow radiation in the zone b > bo ^ 
1/4 is quickly absorbed and transformed into secondary 
pairs. This has a deceleration effect on the outfiow, as 
the energy per particle is reduced (while the net energy 
of particles is conserved, since almost no photons escape 
this zone). 

2.3.2. Radiative zone 

As the outfiow approaches bo, the scattered photons 
begin to escape. In the region 6 < 0.1 all scattered pho- 
tons avoid conversion and splitting, and escape the mag- 
netosphere, forming the hard X-ray component of the 
magnetar spectrum. As the outflow reaches the region 
b ~ 0.01, 7sc decreases to a mildly relativistic value ~ 1 
(see Equation [5]) , i.e. essentially all its energy has been 
radiated away. 

We conclude that the relativistic outflow radiates al- 
most all its energy in a well-defined region of magnetic 
field 0.01 ^ b ^ 1/4, regardless of the details of magneto- 
spheric configuration — the same conclusion is valid for a 
multipolar twisted magnetic field. Section 3 will describe 
a simple method to calculate the produced radiation. 

2.4. Multiplicity of pairs 

Let 2N be the number fiux of outfiowing along both 
legs of the magnetic loop toward its top. Multiplicity is 
defined by 

eiV 

X = ^, (9) 

where / is the electric current circulating in the magnetic 
loop. N grows in the pair-loading zone 6 > 6o to its final 
value that can be evaluated at 6 « 6o using 
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Here 7(&o) is the average outfiow Lorentz factor at 6 w 60. 
It is close to 7sc given in Equation ([5]) ; 7 ^ 7sc is expected 
if the particle momentum distribution has a significant 
width (see Section 2.5). Then one finds. 
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where we used the estimate 7(&o) ^ 7sc(^o) ^ 20. 

The large multiplicity of created pairs implies that the 
outfiow can easily conduct the electric current / required 
by the twisted magnetic field and screen electric fields. 
Therefore, no "gaps" with accelerating electric field are 
expected in the outer magnetosphere. 

2.5. Momentum distribution and collective effects 

A simplest model of outfiow would assume that all 
particles have the Lorentz factor given by Equation ([5]). 
In fact, the local distribution of Lorentz factor cannot be 



a delta- function i5(7 — 7sc) — the real distribution must 
be broad. It is infiuenced by three effects: 

(1) Cold, single-fiuid fiow would be unable to carry the 
required electric current. An electric field is induced in 
the magnetosphere to separate the average velocities of 
the positive and negative charges, so that the electric 
current is organized in the (nearly neutral) outflowing 
plasma (B12). It enforces a minimum width of the mo- 
mentum distribution of e^. 

(2) The cascade in the pair-loading zone injects particles 
with a broad range of momenta (see Section 3.2 below). 

(3) The momentum distribution of created particles fe{p) 
is prone to two-stream instability. The instability tends 
to thermalize the distribution and reduce the gradients 
dfe/dp that are responsible for the instability. 

Remarkably, it is possible to evaluate the outflow radi- 
ation even without knowing exactly how /g (p) is changed 
by collective electric flelds. This is possible because 
the outfiow with a broad distribution function still must 
scatter radiation with Lorentz factor 7sc given by Equa- 
tion (O, and the scattering rate is still dictated by how 
quickly the outfiow must decelerate, keeping the parti- 
cles at 7 7sc and avoiding excessive drag. The surfing 
described in the end of Section 2.2 is still a valid descrip- 
tion of the outfiow dynamics. A large width of fe{p) only 
implies that at any given location a fraction of particles 
actually scatter radiation — those particles that have 
7 « 7s(. — while all other particles do not have suffi- 
ciently high Lorentz factors for resonant scattering. The 
energies of scattered photons Esc and the net scattering 
rate are not much different from the single-fluid model. 
This fact is observed in the simulations of B12, where 
fe{p) is broad. 

3. CALCULATION OF HARD X-RAY EMISSION 

We build the emission model starting from the sim- 
ple consideration of a single particle. Then we consider 
emission from the cascade that develops along the 
magnetic loop, without taking into account collective ef- 
fects. Finally, we formulate an approximate model for 
the outflow with collective behavior, which will be 
used to calculate sample spectra in Section 4. 

3.1. Emission from one particle 

Consider the simplest case of one particle shown in Fig- 
ure 2. The particle is injected with a high Lorentz factor 
in the low part of the magnetic loop and decelerates as 
it moves away from the star. What is the net spectrum 
of radiation emitted by the particle as it decelerates to 
7~ 1? 

The Lorentz factor 7 decreases along the "scattering 
curve," 7 ~ 7sc, which is given by Equation ([5]). It will 
be convenient to rewrite Equation ^ in the form. 



-fsc{b) = Kb, 
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The particle energy is passed to high-energy photons that 
are generated by scattering thermal radiation. Each scat- 
tered photon has energy E = brUeC^ in the particle rest 
frame, as the scattering is resonant^ Its blueshifted en- 
ergy in the observer frame E depends on the scattering 

^ We consider the radiative zone where 6 <g 1 and therefore 
neglect corrections due to electron recoil. 
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angle. A simplest estimate for the emitted spectrum is 
obtained if we neglect the dispersion in scattering angles 
and take E equal to its average value jE, 

E^jbmeC^. (13) 

Since the particle moves with 7 Kb (Figure 2), we 
can substitute b = j/K in Equation (jl3p and obtain the 
relation between the emitted photon energy E and the 
particle energy jrrieC^, 

7 ruec^ = {K rrieC^E) . (14) 

Energy dWi = dj nieC^ is radiated in photons of energy 
E, and the net emitted spectrum is given by 

2\ 1/2 

■ (15) 
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Since E (xb^, the most energetic photons are emitted at 
the highest b ^ bo, near the boundary of the radiative 
zone, with energies comparable to 1 MeV. As the particle 
moves to smaller magnetic fields b and its Lorentz fac- 
tor decreases, softer photons are produced, forming the 
spectrum oc E^^^^ . 

3.2. Emission from cascade (no collective effects) 

Again consider one high-energy particle injected with 
a Lorentz factor 7p > 100 at 6 > 1. Emission from the 
particle in the radiative zone b < bo was evaluated in 
Section 3.1. Before emitting the escaping radiation, the 
particle scatters photons in the zone b > bo and thus gen- 
erates secondary e"^ pairs. The secondary particles also 
scatter photons, which can convert to more pairs. All 
pairs generated in this cascade enter the radiative zone 
and eventually decelerate to 7 ^ 1 as they outflow along 
the magnetic loop. What emission should be observed 
from the cascade? 

In this section, we adopt the assumption that all parti- 
cles flow freely and change their 7 only when they scatter 
a photon. This model ignores the interaction between 
particles via collective electric fields (Section 2.5). The 
opposite case of strong collective effects will be consid- 
ered in Section 3.3. 

The development of the cascade is shown in Figure 3; 
the details of the calculation are given in Appendices A 
and B. Each created particle follows a track in the B-j 
plane similar to that in Figure 2, except that the sec- 
ondary particles are injected with smaller energies. They 
start their tracks significantly below the scattering curve 
7sc = Kb and have a long free path with 7 — const before 
reaching the curve; then the particles begin to frequently 
scatter radiation and surf along the curve as shown in 
Section 2.2. 

The cascade develops in the zone of a strong mag- 
netic field where the photon splitting rate is high (Ap- 
pendix B). The scattered photon can be in one of the 
two polarization states, _L or |j. The energetic _L photon 
quickly splits in two daughter photons with || polariza- 
tion, before it gets a chance to convert to e^. The | 
photons do not split; they convert to e^. The conversion 
occurs immediately after the propagating photon meets 
the threshold condition ([5]). Therefore, and e+ are 
created with total energy E± w i?thr, in the ground Lan- 
dau state, and do not produce synchrotron emission. 
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Figure 3. One random realization of cascade produced by 
one electron injected with Lorentz factor 7p = 10^ at magnetic field 
B = WBq . 136 secondary particles are created in this cascade; the 
injection points for the primary particle and 68 secondary pairs are 
shown by blue crosses. After injection, each particle moves toward 
and then along the "scattering curve," which is approximately de- 
scribed by Equation (0. 

Note that pair creation takes place in the region 6 > 0.1 
where photons should convert to positronium rather than 
free (Shabad & Usov 1986). The binding energy of 
positronium atoms is Et ^ 0.1 keV, and they are eas- 
ily unbound by photons of energy comparable to Eb (in 
the atom rest frame), i.e. by UV photons of energy 
^ 0.l7~^ keV in the lab frame, where 7 = E±/2meC^ 
is the Lorentz factor of the positronium atom. The rate 
of this reaction in a thermal radiation field was evaluated 
by Bhatia et al. (1988). In our problem, the positron- 
ium Lorentz factors are modest, and the resulting ioniza- 
tion rate is so high that essentially all pairs must be 
unbound. Even if a positronium atom reaches the scat- 
tering curve before ionization, and e~ begin to reso- 
nantly interact with X-rays as if they were free particles 
(because huB 3> Et) and immediately become unbound 
due to recoil in scattering. 

Our cascade simulation shows that the total number of 
secondary pairs (generated per one electron injected 
with Lorentz factor •jp = 10'^) is TVinj ~ 10^, similar 
to estimate (fTT|) . The particle Lorentz factors at injec- 
tion points (blue crosses in Figure 3) form a distribution 
dNinj/d'j that can be accurately calculated by simulating 
a large number Af of random realizations of the cascade. 
We use J\f — 10* to accumulate sufficient statistics. The 
obtained distribution dA'inj/c?ln7 is shown in Figure 4. 
It is very well approximated by a simple power law. 
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7 > 71. 
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Here 71 = E/2meC^ ~ 1 — 2 corresponds to photons of 
minimum energy E that are able to convert to e^; in our 
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Figure 4. Injection Lorentz factor distribution of pairs in 
cascades generated by one particle with 7p = 10'^ (of. blue crosses 
in Figure 3, which shows one realization of the cascade). The 
injection of the primary particle is shown by the peak at 7 = 10'^ 
(normalized to unity), and the broad smooth distribution shows the 
injection of second ary e^. Straight line shows the approximation 
given in Equation I II6II . 



sample model 71 = 1.5. The slope of the distribution in 
Figure 4 is g w 0.84, and the corresponding normaliza- 



tion factor A is given by ^ = q^lN^^ 
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A cascade is characterized by its "pair yield" Y . It 
is defined as the mean expectation for the ratio of the 
total rest-mass energy of created pairs, A'injTOeC^, to the 



energy of the primary particle 7^7 



We find 



Y 



0.15. 



(17) 



The same Monte-Carlo simulation gives the spectrum 
of radiation emitted by the cascade. We accumulate the 
statistics of all escaping photons (which avoid splitting 
and conversion) in A/" = 10"* realizations of the cascade, 
and find their average spectrum. The result is shown in 
Figure 5. In this simulation, we assume that radiation is 
produced by a thin magnetic tube for which the highest 



energy of escaping photons E-^ 



SrrieC (these pho- 



tons are emitted at & w &o ~ 0.3); therefore the obtained 
spectrum has a sharp break at £^ fa 1.5 MeV. We also 
performed more detailed simulations that track all pho- 
tons as they propagate through the magnetosphere and 
check for absorption or splitting along their trajectories 
(Appendix B); this gave a similar break at -Emax whose 
exact value (comparable to SrUeC^) depends on the radius 
of curvature of the magnetic tube. 

A simple analytical estimate for the spectral slope be- 
low 1 MeV can be derived as follows. Each injected parti- 
cle begins to radiate when it reaches the scattering curve 
7 = Kb. Each particle surfing along the curve emits 
the same standard spectrum dWi/dE given by Equa- 
tion (fT5|) . Hence radiation produced at a given photon 
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Figure 5. Average spectrum of escaping radiation from cascades 
generated by primary particles with iLorentz factor 7p = 10^. The 
total emitted energy W equals the injected e nerg y ■ypnieC^. Dashed 
line shows the analytical estimate (Equation I20I I. Photon splitting 
creates an additional bump that is visible at iJ ~ 1 MeV. The 
bump consists of photons with || polarization. Emission at lower 
energies i? <C 1 MeV is dominated by photons with ± polarization. 



energy E is given by 

dW _ dWi dNinj 
llE " ~dF d-f 



KE 



1/2 



where 7^ is the Lorentz factor of particles that emit at 
energy E (Equation [T^. From Equations (fTS)) and ([11" 
we find 



dW A 



d\nE 2q 
which gives 



— nieC 



dW 
d\nE 



KE 



OAjpirieC^ 



E 



ITleC 



2 > I' (19) 



KE 



0.1 



(20) 



where we used Equation (fT7|) . A/q = A^inj7i , ^^'^ 9 ~ 0.8. 
This expression only describes radiation produced by 
scattering; it does not take into account photon splitting. 
It gives a reasonable approximation to the emitted spec- 
trum at < nieC^ ■ The splitting of high-energy photons 
with the _L polarization, -L— )-|| -I- ||, creates additional 
photons that form a bump near E ^ 1 MeV (Figure 5). 

The cascade model described above has one feature 
that makes the model questionable: the momentum dis- 
tribution of develops two peaks and becomes unstable. 
We have calculated the distribution of created by the 
cascade and how it evolves with b as the flow moves along 
the tube. The momentum distribution at three locations 
in the tube is shown in Figure 6. One can see the devel- 
opment of two peaks. The low-energy peak (p « 1.5 at 
b = 0.07) is formed by recently injected secondary parti- 
cles. The high-energy peak (p 7 at 6 = 0.07) is formed 
by older particles that were injected at higher b and are 
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Figure 6. Momentum distribution of outflowing pairs cre- 
ated by the cascade, without taking into account collective effects. 
The distribution changes along the magnetic tube with decreas- 
ing b = B/Bq as more particles are created and the particles lose 
momentum to scattering. 

now moving along the scattering curve (cf. Figure 3); 
its position evolves with b according to p « ~ 1006 
(Equation [T2|) . 

The two-peak distribution is prone to plasma insta- 
bility that excites Langmuir waves in the outflow. We 
conclude that the model without collective effects is un- 
stable, and the excitation of collective electric fields is 
inevitable, which will change the distribution function. 
Below we will seek a simplest emission model for the 
outflow with collective behavior. 

3.3. Emission from outflow with collective behavior 

At a given location in the active loop, only a small 
fraction of particles have 7 ~ jsc that is high enough to 
resonantly scatter radiation; most particles have 7 < 7sc 
and do not interact with radiation. This fact is true in 
the absence of collective effects (see Figure 3) and re- 
mains true in collective models (cf. numerical simula- 
tions in B12). Since radiation is scattered by particles 
with known Lorentz factor 7 « 7sc (Equation I12|) , the 
details of /e {p) are not important for the produced emis- 
sion. It turns out that it is sufficient to know the average 
Lorentz factor 7(&), 



7= / {i+p'y^'up)dp. 



(21) 



Our model of outflow with collective behavior will 
avoid the complicated calculation of /e(p); instead, we 
focus on 7 and use the following approximation, 

7 0C7sc, (22) 

where 7sc is given by Equation p^ . In essence, j{b) is 
kept somewhat below jsc{b) — the function jsc{b) shows 
how the outflow must decelerate as it enters weaker mag- 
netic fields. As an example of collective behavior, one can 



picture an outfiow whose distribution function maintains 
a constant shape, which shifts to lower 7 as the outfiow 
moves to weaker magnetic fields. Approximation (|22p 
greatly simplifies the calculation of generated radiation. 
Remarkably, the details of the cascade, its multiplicity, 
and the coefficient of proportionality in Equation (|22|) 
turn out unimportant — they drop out from the calcu- 
lation as will be seen below. 

Consider a thin active magnetic fiux tube. We wish to 
know how its luminosity L is distributed over b along the 
tube, as this will determine the spectrum of produced 
radiation. We will use the simple fact that the power 
emitted by the tube element of length dl equals the power 
lost by the outflow, i.e. the emitted power is determined 
by the deceleration law. 

In the radiative zone, where no new pairs are cre- 
ated, the particle number flux TV is constant along the 
tube up to its summit, where particles annihilate. The 
outfiow power 2N^me<? (factor of 2 for two legs of the 
tube) decreases proportionally to 7. The lost power 
transforms to radiation, and hence the produced lumi- 
nosity is distributed over b according to 



dL „ ,v d'y , 

— = 2iV— f TOeC = const. 

db db 



(23) 



Here b — B/Bq is used as a coordinate along the mag- 
netic fiux tube; expression on the left side is constant 
along the tube because 7 is proportional to 7sc (Equa- 
tion and 7sc is proportional to b (Equation [T^ , so 
d'y/db — const. The total power L (emitted in the radia- 
tive zone b < bo) equals {dL/db)bQ, which gives 



dL 



L 
b'o' 



(24) 



Equations ([T^ and ([^ lead to a complete description 
of the emitted radiation, as at every b they determine the 
Lorentz factor of the emitting particles 7sc and the emis- 
sion rate. Consider first the net spectrum of radiation 
emitted by the tube in all directions. Similarly to Sec- 
tion 2.1, a simplest estimate is obtained if we neglect the 
dispersion in scattering angles and take photon energy E 
equal to its mean expectation jscbrrigC^, 



E = Kmy b^ 



(25) 



Then one finds db {l/i){K nieC^ E)-^/'^ dE, and sub- 
stitution to Equation (IMl) gives the spectrum emitted by 
the tube, 

dL _ L 

dE ~ 2bo (KnieC-^Ey/^' 

The total produced luminosity is proportional to 
J E^^^^ dE and peaks at the high-energy end as E^^^. 
The most energetic escaping photons are emitted at 
b ~ bo, near the boundary of the radiative zone, with 
energies comparable to 1 MeV. Equation p6| does not 
include the additional component that is generated near 
1 MeV by photon splitting; it will be included in the 
detailed Monte-Carlo models presented below. 

In Section 2.3, the cascade without collective effects 
gave a softer spectrum dL/dE cx £;^0-9 (Equation [201) • 
In that case, 7(6) was not proportional to 7sc; the sec- 
ondary particles with 7 < 7sc kept constant 7 until they 
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reached the scattering curve at lower b and radiated their 
energy in softer photons. By contrast, the collective be- 
havior with 7 oc invokes a quick redistribution of 
energy between the injected particles with 7 < 7sc and 
the radiating particles with 7 « ■ 

3.4. Detailed calculation: angular distribution 

To evaluate the emission received from an active mag- 
netosphere by observers with various lines of sight, we 
need a more detailed calculation. A concrete geometry 
of the j-bundle must be chosen in such simulations. 

One can view the active part of the magnetosphere 
(j-bundle) as a collection of infinitesimally thin active 
magnetic flux tubes. The flux tubes are assumed to ex- 
tend sufhcicntly far from the neutron star, reaching the 
region oi B < 10^^ G. Consider one elementary tube 
that contains magnetic flux Sr and electric current 61. 
The discharge operates in the tube under voltage $ and 
creates an outflow with power. 



SL = x^SI. 



(27) 



Here we allowed a fraction 1 — x of the total ohmic power 
released in the discharge zone to go to the footpoints of 
the tube, and the fraction x is given to the outflowing 
particles. 

Emission from each piece of the tube is beamed along 
the outflow, and the emerging radiation is strongly 
anisotropic. Consider a short piece of the magnetic tube 
of length dl which corresponds to db = (db/dl) dl. Its 
contribution to the tube luminosity is {SL/2bo) \db\ (see 
Equation [M] and take into account that there are two 
elements dl that correspond to db, as the tube has two 
footpoints). The average energy of emitted photons is 
jscbrUeC^ , and hence the photon emission rate for the 
tube element dl is given by 



diSNpy,) 



6L 



\db\ 



(28) 



27sc6TOeC^ bo 

It remains to determine the spectrum and polarization 
of the emitted photons. In the radiative zone b is small 
(6^1) and hence electron recoil per scattering is small. 
In this regime, the description of polarization and angu- 
lar distribution of scattered photons is simple (see Ap- 
pendix A in BI2 for details). Photon energy in the ob- 
server frame, E, is related to its emission angle in the 
electron rest frame, z?, by the Doppler transformation of 
E = Hlob, 

E^hwB7scil + /3scCosd). (29) 

Approximately 3/4 of s catte red photons have the _L po- 
larization (cf . Equation IA3|) ; these photons are emitted 
with an isotropic distribution. 



P± (cos 1)) — const — - . 



(30) 



The remaining 1/4 of photons are emitted with the 
polarization, and their angular distribution is 



(31) 



Using the Monte-Carlo technique, we generate the angu- 
lar distributions for photons with _L and || polarizations. 



and determine the photon energies according to Equa- 
tion ((29t . 

This straightforward Monte-Carlo simulation gives the 
spectrum and angular distribution of emission from all 
parts of the magnetic tube. We determine the fate of 
emitted photons by tracking their trajectories through 
the magnetosphere and checking for splitting or conver- 
sion to , which allows us to smoothly describe the tran- 
sition from pair- loading to radiative zone (Appendix B). 
Photons that avoid splitting and conversion form the ob- 
served spectrum. Finally, the net observed emission is 
obtained by integrating the emission along the tube and 
then taking the sum of emissions from all tubes in the 
j-bundle. 

The same calculation may be performed by integrat- 
ing emission over the volume of the j-bundle in usual 
spherical coordinates. This method is described in Ap- 
pendix C, where emissivity per unit volume dL/dV is 
derived. Appendix C also describes the concrete case of 
a twisted dipole magnetosphere. 

3.5. Multiple scattering of hard X-rays 

The Monte-Carlo simulation described above neglects 
the fact that the generated high-energy photons can scat- 
ter multiple times in the outflow before escaping. It 
turns out that this additional scattering weakly affects 
the emerging spectrum. Note two facts: (1) Hard X-rays 
escape the scattering region on the light-crossing time re- 
gardless of the number of scatterings; this is because the 
plasma moves relativistically and the scattered photons 
are beamed along the flow, so they keep moving away 
from the star. (2) The outflow significantly loses energy 
only when it scatters soft target photons that propagate 
at large angles with respect to the outflow direction; 
then scattering boosts the photon energy by the large 
factor of ~ 7f(. 10 Since the scattering reduces "i^ to 
^ (which corresponds to a quasi-isotropic distribu- 
tion of photons in the outflow rest frame), subsequent 
scatterings have a small effect on the photon energy. 

To check the effect of multiple scattering on the hard 
X-ray spectrum we used the Monte-Carlo transfer code 
developed in the accompanying paper (B12). We fol- 
lowed the propagation and scattering of the generated 
hard X-rays through the flowing magnetospheric plasma 
and found that multiple scattering changes the emerging 
spectrum by less than 20%. This change makes the spec- 
trum slightly harder; it will be neglected in this paper. 

4. SAMPLE NUMERICAL MODELS 
4.1. Model setup 

Consider a neutron star of radius i? = 10 km, tem- 
perature kT = 0.5 keV, and magnetic dipole moment 
fi = 5 X 10^^ G (which corresponds to the magnetic 
field at the pole -Bpoio = 10^'' G). Poloidal magnetic field 
is well approximated as dipole if the twist amplitude is 
moderate, ^ ^ 1 (B09; Parfrey et al. 2012b). We will use 
spherical coordinates r, 9, <j) with the polar axis aligned 
with the dipole moment; = on the axis. 

Consider an activated magnetic loop (j-bundle) whose 

This boost happens in the lab frame; the photon energy is 
unchanged in the outflow rest frame — only its angle is changed 
by scattering (from cos 1? ~ — 1 to a random value) . 
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northern footpoint on the star is defined by 

y < < 0., < < 0.. (32) 

In our sample models, sin^ 6*^ 0.1; is left as a free 
parameter (0* = 27r describes an axisymmetric j-bundle) . 
Practically the same results are obtained if the j-bundle 
includes the entire polar cap < 6* < 6**, as the field 
lines with footpoints a.t 9 < 0^/2 contribute only a small 
fraction ^ 2^^ to the emitted luminosityO The active 
loop extends to radii r ^ R/ sin^ 6^, lOR. We wish to 
know the spectrum of radiation emitted by the loop, as 
viewed along any given line of sight. 

The total luminosity of the loop is controlled by its 
twist amplitude -0 and the discharge voltage $ (Sec- 
tion 2.1). It is approximately given by (see B09 and 
Appendix C) 

c$/i?A sin'' 6** 
^ AR^ [2^)- ^^^^ 

The discharge occurs in both hemispheres, forming two 
outflows that meet and annihilate at the equatorial plane 
of the magnetic dipole. Therefore, the total emitted 
power corresponds to $ that is two times larger than 
the threshold for the discharge; $ = 10^ — lO^*' V is ex- 
pected (Beloborodov & Thompson 2007) and consistent 
with observations (Beloborodov 2011). It is convenient 
to define the parameter 

Then the total power emitted by the loop may be written 
as 

The observed emission from the loop is calculated using 
the Monte-Carlo technique as described in Section 3.4. 
We assume that the outflow behaves collectively (Sec- 
tion 3.3), as this model is both simpler and more realistic 
than the cascade without collective effects. 

The emission from the loop is strongly anisotropic and 
hence emission observed along a given line of sight will be 
modulated by the rotation of the object. Let us choose 
the fixed (non-rotating) lab frame so that its z-axis is 
parallel to the spin axis, and the xz plane contains the 
observer's line of sight. The instantaneous orientation 
of the magnetosphere with respect to the lab frame is 
described by three Euler angles a, ip, C- Here a is the 
angle between the magnetic and rotation axes, if is the 
angle of rotation about the spin axis, and — ^ = Anodes is 
the azimuthal angle (measured in magnetic coordinates 

r, 9, (j)) of the line of nodes Q x fi. 

Let /3 be the angle between the line of sight and the 
rotation axis. The instantaneous observed emission is 
determined by four angles a, ip, C, and /?, and the ob- 
served spectral luminosity Le varies periodically with ip. 
The magnitude of the variation may be described by the 

* Electric current in the j-bundle 6 < 6t with a given twist i/i 
scales as / oc i/i sin* 6^, (Appendix C). 



"pulsed fraction," 

T max T min 

UE) = ^ (36) 

E E 

or the "pulsed spectrum," 

7- puis 7- max 7- min /'Q'7^ 

^E = ^E ^ ^E ■ l-J'J 

If the emission is not resolved in rotational phase, its 
observed spectrum is obtained by averaging the instan- 
taneous Le over the rotation period, 

1 /-^^ 

Le = — Le dep. (38) 
27r Jo 

The result depends on a, /3, and (. If the emitting re- 
gion (j-bundle) is symmetric about the magnetic axis, 
the result only depends on a and /3. 

4.2. Aligned rotator 

Aligned rotator has a = (magnetic axis parallel to 
the rotation axis) . Then (p = 4> and averaging over rota- 
tion is the same as averaging over the magnetic azimuthal 
angle The shape of the averaged spectrum Le does 
not depend on the azimuthal extension of the j-bundle 
(pi,; e.g. one may chose (p^, — 2-k (axisymmetric j-bundle). 

The spectrum produced by the axisymmetric j-bundle 
is shown in Figure 7. It strongly depends on the angle 
between the line of sight and the magnetic axis. In par- 
ticular, the position of the high-energy peak significantly 
varies around 1 MeV. One can qualitatively understand 
the variations in the spectrum with changing inclination 
by combining the following facts. (1) Emission from the 
loop is generated on field lines that have footpoints at co- 
latitudes 6 ^ 0.3 and reach the equatorial plane at radii 
^max — 10—20. (2) Emission is strongly beamed outward 
along the magnetic field lines. (3) Photons of a given en- 
ergy E are mainly generated where b « 0.1(E /mgC^)^^^ 
(see Equation [^5)1 . Note also that observers in the equa- 
torial plane (inclination 9 — 90°) are special — by sym- 
metry, both hemispheres equally contribute to the ob- 
served spectrum. At small inclinations, one hemisphere 
strongly dominates hard X-ray emission; however, the 
contribution of the other hemisphere is still significant in 
the soft-X-ray band. 

If the emission shown in Figure 7 is averaged over co- 
sine of inclination angle one obtains the spectrum shown 
in Figure 8. Below 1 MeV, it approximately follows the 
power-law ELe oc in agreement with Equation (|26p 
derived in Section 3.4. It corresponds to photon index 
F = 1.5 (which is defined by Le oc E-'^+'^). This does 
not imply that, given a random inclination, the most 
probable observed spectrum has F — 1.5. Significantly 
harder slopes F ~ 1 are observed in Figure 7. In a range 
of inclinations, the spectra have breaks well below 1 MeV, 
and after averaging over inclination this gives the feature- 
less power law with a softer slope F = 1.5. 

Approximately 3/4 of photons scattered in the j-bundle 
have the _L polarization. Therefore at energies E < 
1 MeV (where the splitting of _L photons is negligible), 
the produced emission is dominated by the _L polariza- 
tion, see the lower panel in Figure 8. 

The inclination-averaged spectrum exhibits a signifi- 
cant hump at i? ^ 1 MeV. It is the result of photon split- 
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Figure 7. Spectrum from the axisymmetric j-bundle, observed at 
four different angles with respect to the magnetic axis: 15°, 30°, 
60°, and 90°. The spectrum is norm alize d to the total luminosity 
of the j-bundle L given by Equation J35}. 
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Figure 8. Upper panel: spectrum from the axisymmetric j- 
bundle, Lg, averaged over the object inclination to the line of 
sight. Lower panel: relative contributions of -L and || to the total 

spectrum = Lj^ + 

ting. The daughter photons produced by splitting have 
the 11 polarization, and the splitting-dominated compo- 
nent is particularly visible in the lower panel of Figure 8 
that shows the polarization of emitted photons. Note 
that Figure 8 does not take into account the possible 
multiple scattering of high-energy photons on their way 
out of the j-bundle (Section 3.5). Additional scattering 
weakly changes the spectrum of escaping radiation, how- 



ever it can significantly affect polarization of the MeV 
hump. Scattering of daughter photons after splitting 
switches polarization from || to _L with probability of 3/4 
and tends to bring the polarization to the same value as 
below 1 MeV. 

No pulsations are produced by the aligned rotator if 
the j-bundle is axially symmetric, (f)i, = 27r. Strong pul- 
sations are produced if < 27r, i.e. if the j-bundle is 
confined to a smaller range of azimuthal angles. The pul- 
sations are the consequence of strong beaming of high- 
energy emission along the magnetic field lines, which 
makes the emission almost invisible when the line of sight 
is outside the interval < < 0*. The pulse width is 
proportional to (/i^, and the pulse fraction can approach 
100%. 

4.3. Orthogonal rotator 

Orthogonal rotator has a = t: /2 (magnetic axis is per- 
pendicular to the rotation axis). Then even an axisym- 
metric j-bundle (0* = 27r) will produce pulsating emis- 
sion. As the magnetar rotates, the inclination of line of 
sight to the magnetic axis, 0, periodically changes and 
the observer samples the spectra shown in Figure 7. The 
spectrum strongly depends on 9 and hence significant pe- 
riodic variations in are expected, depending on /3 (the 
angle between the line of sight and the rotation axis). In 
particular, the position of the high-energy peak will os- 
cillate with a significant amplitude. No pulsations are 
expected only if /? = 0; the observed spectrum in this 
case is the same as the 9 = 7r/2 spectrum in Figure 7. 

Figure 9 shows the emission from the orthogonal ro- 
tator viewed at /9 = 20°. The rotation-averaged spec- 
trum has the photon index F « 1 at photon energies be- 
low 70 keV. The pulsed spectrum is significantly harder, 
F « 0. The pulsed fraction increases with E and ap- 
proaches 100% aX E K, 300 keV. The model is close to 
observational data reported for 4U 0142-1-61 and IRXS 
J170849-400910 — two magnetars whose hard X-ray pul- 
sations have been studied in detail (den Hartog et al. 
2008a,b). 

Figure 10 shows the emission from the orthogonal ro- 
tator viewed at /3 = 90°. The rotation-averaged spec- 
trum has the photon index F « 1.5 below 1 MeV; this 
spectrum is close to the average spectrum in Figure 8, 
although not exactly the samejj The pulsed spectrum is 
again harder than the average spectrum. The amplitude 
of pulsations is significantly larger than in Figure 9; the 
pulsed fraction approaches 100% at i? « 20 keV. 

The variety of observed spectra and pulse profiles be- 
comes even larger if the j-bundle is not axisymmetric. 
Then the observed emission depends on the third Euler 
angle C (Section 4.1). The pulse profile in a given spectral 
window depends on the angles a, /?, C f^nd the geometry 
of the j-bundle. 

5. DISCUSSION 

In this paper, we calculated the X-ray emission ex- 
pected from large twisted magnetic loops around mag- 
netars. The loops are inevitably heavily loaded with 

^ For the orthogonal rotator viewed at /3 = 90°, we have d = ip, 
where d is the angle between the line of sight and the magnetic axis. 
Therefore, averaging over rotation < < 27r is equivalent to 
averaging over d. In Figure 8, the emitted spectrum was averaged 
over cosd (rather than d), giving a slightly different result. 
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Figure 9. Spectrum from the orthogonal rotator (a = 90°) when 
the line of sight is at angle /3 = 20° with respect to the spin axis. 
The j-bundle is assumed to be symmetric about the magnetic axis 
(0* = 27r). So lid curve shows the spectrum averaged over rotation 
(Eq uati on I38I I. Dashed curve shows the pulsed spectrum (Equa- 
tion |33. 
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Figure 10. Same as Figure 9 but for /3 = 90°. 

pairs, with multiplicity ~ 10^, which leads to efScient 
screening of electric fields. The created plasma outflows 
along the magnetic loop and loses its energy to radiation. 
Transformation of the outflow energy to escaping radia- 
tion occurs in the region lO^'^ > i? > 3 x 10^^ G, which 
requires a minimum size of the loop 3 to 10 stellar radii 
(if the surface field is ~ 10^^ G). We did not consider 
small loops because most of photons scattered near the 
star are unable to escape the magnetosphere; besides, 
smaller loops are less likely to be active (Section 2.1). 

A key feature of the outflow along the loop is that 
its Lorentz factor is self-regulated so that it scatters ra- 



diation with 7sc w IQQB/Bq (Section 2.2). This fact de- 
termines the radiation spectrum emitted by the outflow, 
taking into account that the scattered photons have the 
characteristic energy E ~ ^sc{B / BQ)meC^ . The emit- 
ted spectrum below 1 MeV is hard, i.e. softer photons 
E <^ \ MeV contribute less to the outflow luminosity. 
This is simply because softer photons are emitted where 
the outflow decelerates, and there remains less energy 
to emit. The spectrum is also suppressed at energies 
£^ ^ 1 MeV, because photons of such high energies are 
emitted in the zone of strong B and cannot escape. 

The model predictions appear to agree with available 
data. Observations of magnetars in the 10-100 keV band 
show a typical photon index F ~ 1 — 1.5, similar to what 
is found in the model spectra. The observed spectra 
peak above 100 keV, and the existing upper limits for 
AXP 4U 0142+61 at 1-30 MeV indicate a spectral break 
between 0.3 and 1 MeV (den Hartog et al. 2008a). 

The strong beaming of radiation along the loop ex- 
plains why the observed hard X-ray flux shows huge pul- 
sations, with the pulsed fraction growing with photon 
energy up to 100% (Kuiper et al. 2003; den Hartog et 
al. 2008a, b). No unique pulse profile is predicted by our 
model, because the pulse is sensitive to the loop geome- 
try, its orientation with respect to the rotation axis and 
the line of sight. The pulse profile may be further com- 
plicated if more than one loop are activated and the su- 
perposition of their emissions is observed. The expected 
hard X-ray pulse does not have to be in phase with the 
soft X-ray pulse, as observed in AXP IRXS J170849- 
400910 (see Figure 6 in Den Hartog et al. 2008b). 

In general, sources of beamed, pulsed radiation are not 
expected to display a standard spectrum at all viewing 
angles, and magnetars are not an exception. The active 
loop typically emits a hard spectrum, F ~ 1 — 3/2, but 
its detailed shape depends on the loop and the viewing 
angle (Figures 7, 9, 10). The predicted high-energy peak 
significantly depends on the line of sight and its observed 
position is expected to oscillate as the magnetar rotates. 

Note also that a very hard spectrum F < 1 may be 
produced if smaller magnetic loops are activated. For 
example, consider a loop that is confined to the region 
of strong magnetic field B > 3 x 10^^ G, where Hujb > 
35 keV. The loop generates photons of energy E ^ ^Tiljb 
and almost no emission is produced below 35 keV, leading 
to a hard spectral slope. 

The surface transition layer between the corona and 
the star may give a second contribution to the hard 
X-ray emission (Thompson & Beloborodov 2005; Be- 
loborodov & Thompson 2007). This dense layer is bom- 
barded by the particles from the discharge and may 
reach temperatures kT <; 100 keV, sufficient to emit 
high-energy bremsstrahlung radiation with photon in- 
dex F w 1. Observed magnetars show complex spec- 
tral variations with rotational phase, which suggests that 
the emitting plasma occupies an extended region in the 
magnetosphere rather than a heated spot on the surface. 
This favors the outfiow model described in the present 
paper. Polarization measurements would provide a fur- 
ther test discriminating between the two mechanisms. 
The surface transition layer emits photons with the | 
polarization, while the loop itself emits radiation that is 
dominated by the _L polarization (also called "extraor- 
dinary" or E-mode) . Careful calculations of polarization 
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should include transfer effects (Fernandez & Davis 2011). 

5.1. Annihilation line 

The outflowing pairs accumulate and annihilate at 
the top of the magnetic loop with Lorentz factors 7^1 
(Figure 1). This must create an annihilation line at 
E w 0.511 MeV, which we did not show in our sample 
spectra. There is a simple relation between the luminos- 
ity of annihilation radiation iann and the total outflow 
power L (which is equal to the luminosity of the hard 
X-ray component). The outflowing particles enter the 
radiative zone B < 10^^ G with the average Lorentz fac- 
tor 7o ^ 10 — 20 and radiate almost all their kinetic 
energy (70 — l)meC^ in hard X-rays before annihilating. 
Therefore, a fraction of ~ of the total outflow energy 
is converted to annihilation radiation, 



L 



0.1. 



70 



(39) 



The visibility of the annihilation line at 511 keV depends 
on its width. As are decelerated before annihilating, 
a pronounced narrow line is expected. Mildly relativistic 
motions in the annihilation region would make the line 
broad and less visible. Note also that Lann is approxi- 
mately isotropic, while L is strongly anisotropic; there- 
fore the observed ratio Lann/^ will depend on the line of 
sight and can be significantly larger or smaller than the 
mean expectation given by Equation p9p . 



5.2. Temporal behavior of activity 

The nonthermal emission of magnetars must be fed 
by ohmic dissipation of electric currents excited in their 
twisted magnetospheres (Section 2.1). The rate of en- 
ergy release is proportional to voltage established along 
the magnetic field lines, which is regulated by electric 
discharge to $ ~ 10^ V, in agreement with the observed 
luminosity. One expects the untwisting magnetosphere 
and its luminosity to evolve on the ohmic timescale, 
which is comparable to 1 year. Several transient magne- 
tars show evidence for this evolution (Beloborodov 2011); 
XTE J1810-197 is a canonical example. Hard X-ray data 
for transient magnetars are meager. A strong evolution 
of the hard X-ray flux is expected in these objects, which 
may be tested by the upcoming NuSTAR observations. 

Other sources, in particular 4U 0142-1-61 and IRXS 
J170849-400910, remain stable over a long observation 
period (about 10 years). Their hard X-ray emission has 
been studied in detail (den Hartog et al. 2008a, b). The 
different patterns of magnetar activity are not under- 
stood because of our ignorance concerning the pattern 
and speed of the crustal motions that cause the activ- 
ity. The crust is likely to flow plastically and gradually 
twist the magnetosphere. It is unclear whether a quasi- 
steady state should be expected, punctuated by sudden 



outbursts triggered by excessive twisting (Parfrey et al. 
2012b). 

5.3. 1-10 keV emission 

In this paper, we focused on the mechanism for the 
hard X-ray peak in magnetar spectra, and did not ad- 
dress the shape of the soft X-ray (keV) peak. As has 
been known for a long time, the keV peak is different 
from simple blackbody (e.g. Woods & Thompson 2006). 
The 1-10 keV spectrum is explained as a superposition 
of the thermal component and a soft power-law tail with 
photon index F = 2 — 4. The soft tail is distinct from the 
hard X-ray component, although a correlation between 
them was reported (Kaspi & Boydstun 2010). 

Previous models of resonant scattering around mag- 
netars were able to reproduce the soft tail in the 1- 
10 keV band (Fernandez & Thompson 2007; Nobili, Tur- 
oUa, & Zane 2008; Rea et al. 2008). In these mod- 
els, the magnetospheric region where fiws ~ 1 — 10 keV 
{B = 10" - 10^2 Q) is flug(i ^ith a mildly relativis- 
tic plasma (7/? ~ 1). It is usually assumed that the 
plasma is made of counter-streaming positive and neg- 
ative charges, which maintain the electric current with 
particle multiplicity ~ 1; then the scattering optical 
depth can be estimated as t ^ 'ip / (3 ^ 1 (Thompson et al. 
2002). The results of the present paper show that these 
assumptions are problematic. The typical pair multiplic- 
ity M must be comparable to 10^, and charges of both 
signs should stream away from the star, toward their an- 
nihilation at the top of the magnetic loop. The model of 
the soft X-ray tail due to resonant scattering should be 
revised accordingly. 

The high multiplicity alleviates the following problem. 
The particle flux along the magnetic field lines is iV = 
M.I/e where I/e ~ 5 x 10^^ V s""'^ is a typical electric 
current circulating through the region hujB ~ 3 keV (cf. 
Equation [1]) . The observed luminosity in the soft tail 
Ltaii is comparable to lO'^^ erg s~^, and hence the energy 
emitted per particle is given by 



ail 



eL 



tail 



2A^ 21 M 



mow 



L 



tail 



J^Q35 g]-g g 1 



MeV. (40) 



The mildly relativistic particles that are supposed to pro- 
duce Ltaii have energies ~ 1 MeV. Hence the picture of 
counter-streaming charges with ~ 1 requires that the 
emitted energy is 2-3 orders of magnitude larger than the 
particle energy^ In contrast, the model with multiplicity 
M <; 10^ implies that a few MeV per particle is sufficient 
to feed the observed Ltaii and modify the thermal peak. 



This work was supported by NASA NNX10AI72G. 



® One could assume that a voltage ^ 1 GeV is applied to this 
region and responsible for the enhanced emission rate per particle. 
In this picture, radiative drag impedes the electrostatic acceleration 
and keeps the particles mildly relativistic. This "radiatively locked 



flow" is considered in more detail in the accompanying paper (B12). 
Note also that the balance between the radiative and electrostatic 
forces is unstable when the particle interacts with photons below 
the Wien peak of thermal radiation. 
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APPENDIX 
RESONANT SCATTERING 
Cross section, photon energy, and polarization 

The polarization states of X-rays in the magnetosphere are determined by the quantum electrodynamic effect of 
vacuum polarization in a strong magnetic field. It defines two eigen modes for an electromagnetic wave with a 
wavevector k: with electric field oscillating along B x k (_L polarization, sometimes called "extraordinary" or E-mode) 
and with electric field oscillating along k x (B x k) (|| polarization, sometimes called "ordinary" or 0-mode). 

In our problem, the outfiow is ultra-relativistic in the region of strong B, and target photons are strongly aberrated 
to become nearly parallel to B in the electron rest frame. Then the scattering cross section is the same for || and _L 
photons. Let Et = tub he the energy of a target photon measured in the electron rest frame. Resonant scattering cross 
section in this frame is given by 

eB 

CT = 27r^reC(5(a; - ws), u)b = , (Al) 

where r^ = /meC^ ■ Equation (jAl|) can be derived classically when B <C Bq (e.g. Canuto et al. 1971; Ventura 1979). 
It also remains valid at any B, as long as target photons propagate nearly along B in the electron rest frame (e.g. 
Harding & Daugherty 1991). 

Consider an electron moving along a magnetic field line with Lorentz factor 7 = (1 — /?^)^^/^. When the electron 
scatters a photon, the photon energy changes from Et to E. Energy conservation gives (7 — 7/) WeC^ = E — Et where 7/ 
is the electron Lorentz factor after scattering. Resonant scattering may be viewed as photon absorption, which excites 
the electron to the first Landau level Eb = (1 + 2fe)^/^TOeC^, followed by photon emission (electron de-excitation). The 
Lorentz factor of the excited electron 71 is found from ^rUeC^ + Et = jiEb- Assuming that de-excitation on average 
does not decelerate or accelerate the electron along B, the mean expectation for 7^ is 71. Then the average energy of 
emitted photons is £' = {jnieC^ + Et){l — rUeC^ / Eb), which yields (using Et <C jniec'^), 

This expression gives a typical energy of scattered photons; the exact E in each scattering event depends on the 
scattering angle. 

The scattered photons can be in one of the two polarization states. Using the accurate relativistic differential cross 
section, we find that the fraction of scattered photons with the _L polarization is conveniently approximated by the 
following formula, 

^ - (A3) 
cr^ + CTii 4 + 6 

Scattering rate 

Let n^{^,ijj) be the local angular distribution of the photon spectral density [cm^'^ s]. The cross section (|A1[) implies 
the following rate of resonant scattering by one electron, 

= J dn J duj n^5{Lo - w^s), (A4) 

where cjres = 7(1~/?'^^) In our problem, the spectrum of target radiation seen by the electron may be approximated 
by a (diluted) Planckian with temperature T and photon density n, 

u:n^^0A2-^nf{n), y = — . (A5) 

Then N^c becomes, 

iVsc « 0.42 n / f{n) y^^^ dn. (A6) 

1 Kl J exp J/rcs - 1 

The resonance condition defines the range of possible ?/rcs, ymin < 2/res < ymax- When j/min ^ 1, only photons in 
the Wien tail of the Planck spectrum can be scattered. Then the integral in Equation ()A6|) strongly peaks where 
J/ros ~ 2/min, which givcs (usiug f[Vt]dVL ^ y~ldy res), 

^^^^ jkT '^•^^'^y--^ ' y™"- 7(1-/3 cos ^?„,ax)fcT' ^"^^^ 

where ymin corresponds to the largest photon angle z9max with respect to {3. Thermal photons reflected by opaque 
plasma near the top of magnetic loops (Figure 1) have the maximum angle, as they stream against the outflow 
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direction. The exponential sensitivity of A^sc to z^max explains why the reflected photons are the main targets for 
scattering even though their total number density n is much smaller than that of photons streaming directly from the 
star. For photons streaming at smaller angles "i?, j/min would increase, leading to an exponential suppression of TVsc- 
For instance, changing j/min = 7 to j/min = 14 reduces Nsc by the factor of 4ex p(— 7) « 4 x 10^'^. 
The density of reflected thermal radiation that should be used in Equation (jA7p is given by 

4^i?2 2.7fcTc ^ ^O.iy ^1035 ergs- V ^0.5 keVy ^100 km/ ' ^ ' 

Here Lth is the thermal luminosity of the star, 2.7fcT ~ 1 keV stands for the average energy of thermal photons, and 
^ is the fraction of Lth that is intercepted by the slow, opaque plasma at the top of magnetic loops. 

In the calculations in this paper, it is convenient to describe the number of thermal photons scattered by an outflowing 
particle as a function of decreasing magnetic field, 

— = ^ A9) 

dB PcdB/dV ^ ' 

where I is length measured along the magnetic field line, and dB/dl (for dipole field lines) is given in Equation (|C8I) : 
in the simulations shown in Figures 2 and 3 we use a simplified dB/dl = —3B/r. 

Drag exerted by the reflected radiation on a relativistic electron 
The electron loses energy to radiation with rate 

4 « -A^sc {E - Et), (AlO) 

where E is the average energy of scattered photons (Equation IA2I) and Et may be neglected compared with E. This 
gives 

m,c3^ a kT \} VT+2b)- ^ ' 

where a = e^ /he = 1/137. The corresponding drag force acting on the relativistic electron is given hy F ~ E,,/ c. 
The expression for E^ may be written in a different, physically more transparent form, 

Ee~-ca,csntE. (A12) 

Here cJres is the effective cross section for resonant scattering, 

CTres- (l-/3cosi9)^reA, (A13) 

where A — 2Trc/ujB is the photon wavelength in the rest frame of the electron. Th e fact or (1 — /3cosi?) (from transfor- 
mation of cross section to the lab frame) is comparable to unity, nt in Equation (|A12|) is the density of photons near 
the resonant energy ftwrcs- Our electron interacts with photons in the Wien tail of the spectrum, and their number 
peaks near the minimum j/ies = hu^cs/kT, 



Uf « ki —— 

dEt 



0.42<i„e-^™n. (A14) 



Substitution of Equations (IA13P and (jA14[) to Equation (IA12P gives an explicit expression for the drag force J- = Ee/c. 
The estimate for T may be summarized in the following form, 

■g(l-/3cosz9max) ^ -10 — gni7y^i„e"''""°, (A15) 



7 rUeC 4 a K 

where nn — n/10^^ cm"^. The recoil correction factor q{h) — 1 when 6^1 (Equation [T]) ■ The approximate estimate 
given in Equation (jA15[) is sufficient in view of the strong exponential dependence on t/min and the large prefactor 
10(r/i?) ni7 y^jjj. The outflowing particles surf with rF/^rUeC^ ~ \d\\iB/d\x\r\ ^ 3 (Section 2.2), which corresponds 
to ymin ^ 7 and the prefactor ~ 10^. 

PHOTON SPLITTING AND CONVERSION TO ELECTRON-POSITRON PAIRS 

Conversion to 

Consider a photon of energy E propagating at an angle ■& with respect to the magnetic field lines. The photon 
can convert an pair if the threshold condition E'thr < -E is satisfied. The threshold energy depends on the photon 
polarization: -Ethr = 2rneC^/ sini9 for || photons (then both and e~ can be created in the ground Landau state) and 
E'thr = 2TOeC^(l -t- 2&)^/^/ sin?? for _L photons (then at least one of the created particles is required to be in an excited 
Landau state). 
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The rate of conversion depends on the local magnetic field b — B/Bq. If 6 ^ 0.05, one can assume practically 
instantaneous conversion once the thershold condition is satisfied; then the e"^ pair is created with the lowest possible 
energy E- + = Ethr- If ^ ^ 0.05, conversion is significantly delayed and generally occurs when E'thr is significantly 
below E. Then are created in high Landau states, and the absorption coefficient for the photon is given by 
Kabs « 4.3 X 10^6 cxp (-SmecVSfe-Esin??) cm^^ (Erber 1966). 

The model described in this paper predicts that practically all photons that convert to do so in the region 6 > 0.1, 
and very close to their emission points. Then the photon energy E is shared by two particles created in the ground 
Landau state with energies E/2. 



Splitting 

Following selection rules (Adler 1971; Usov 2002), we assume that only one splitting channel is allowed: -L^|l + |1 
(Adler [1971] uses the opposite notation for the polarizations states, -Lo||). The absorption coefficient for a _L photon 
of energy E propagating at an angle "& with respect to the magnetic field lines is given by 



Ksp = 0.12 C (6 sin 1?)^ 



E 



cm 



(Bl) 



where C, = 1 \fh and C ^ 46^^ in the limit of 6 ^ 1 (the exact asymptotic value depends on Esuvd, see Baring 
& Harding 1997). Factor can be calculated numerically (Figure 8 in Adler [1971] and Figure 1 in Baring & Harding 
[1997]). Medin & Lai (2010) and our Monte-Carlo simulations use an approximate formula. 



c 



= (£f + 0.05)-^(0.255 + 20)-\ £f = 6^exp 



-0.6 



Esmd 
2meC^ 



(B2) 



When a photon splits, its energy E is shared between two daughter photons propagating in the same direction. The 
probability distribution for the energy a daughter photon E\ is approximately given by p{u) ~ 30u^(l — u)^, where 
u = Ei/E (Baring & Harding 1997). We use this distribution in our Monte-Carlo simulations. 



Implementation of splitting and conversion in Monte- Carlo simulations 

To determine the fate of high-energy photons generated by resonant scattering we follow the photon along its 
trajectory and track the change of its pitch angle i? with respect to the local magnetic field lines. The (initially small) 
angle grows due to the curvature of the field lines while the photon polarization remains unchanged along the trajectory 
(the small difference of the refraction indices for the two polarization states results in locking/ adiabatic tracking of the 
photon polarization). For photons with || polarization, we check if they convert to e^. The threshold for conversion 
_Ethr = Ira^c? I sin I? is reduced as i? increases. When Eti^r/E is reduced below unity, conversion becomes possible, and 
in practice occurs immediately. For photons with ± polarization, wc check if they split or convert to e^. In practice 
conversion is negligible as splitting occurs before the conversion threshold condition is satisfied. 

Consider a photon emitted along the outflow direction at point rem- Let -Emax(rom) be the maximum energy for the 

photon to avoid destruction and escape. More precisely, wc define two maximum energies i?max(i"em) and -£'max(''em) 

for II and _L photons, respectively. For any given magnetospheric configuration, £'Jiiax(rom) and E^^^{rcin) can be 
calculated numerically by following the photon trajectories as described above. 

The fate of an emitted photon in our Monte-Carlo simulations is determined by comparing its energy E with i?max. 
The simulation can be further simplified due to the following observation: if a resonantly scattered photon is destroyed, 
the destruction happens quickly, close to the emission point rem (except for a negligible number of photons scattered 
near the magnetic axis, where the field line curvature is small). This is because the high-energy photons capable of 
conversion or splitting are generated in the region of strong and curved magnetic field & ^ 0.1 (recall that the energy 
of a resonantly scattered photon does not exceed 2'yscbmeC^ ~ 2006^ mgC^). In the Monte-Carlo simulations it is a 
good approximation to assume the immediate conversion of || photons ii E > £Jmax(rein) and the immediate splitting 
of _L photons \i E > E^^^^{y^,^-^). This approximation is used in all sample models shown in Section 4. 

Let _E(roin) be the average energy of resonantly scattered photons. Most of the photons generated in the zone 
E > -Emax are destroyed. As the outfiow moves along a field line, there is a sharp transition to photon escape at 
the boundary tq definc;d by -E'(ro) = ii^max(ro)- Similar results will be obtained if one assumes that all photons 
scattered with E < -E'max(ro) escape (regardless of their emission points rem), and all photons with E > i?max(ro) 
split or convert, depending on their polarization states. Such an approximation is used in the cascade simulations in 

Section 3.2. Instead of the pre-calculated ii^max(rom) and E^g^^{rcm) (which depend on the chosen dipole magnetic 
configuration), we assume a given fixed i?max — -Emax ~ 7sc(&o) rricC^. The value of E-a^g^^ is typically somewhat above 
1 MeV, depending on the curvature of magnetic field lines. Simulations shown in Figures 3-6 assume i?max = 3meC^ 
— a typical value obtained with the more detailed simulations that follow the propagation of photons through the 
magnetosphere. 
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Boundary between the pair-loading and radiative zones 

The location of the boundary between the two zones was estimated in Section 2.3. Approximately, it is where the 
magnetic field b = B/Bq equals 60 ^ 1/4. The exact 60 varies for different field lines of the magnetosphere. It is 
convenient to label the field lines by their maximum extension radius i?max (the apex radius). The boundary location 
also depends on the magnetic moment of the star fi and its temperature T. These dependencies may be formulated 
once we specify a more rigorous definition of the boundary. 

To define the radiative zone, consider a photon scattered at a given point P by an electron with Lorentz factor 7sc 
(Equation [T^ . We focus on the region where b 1 and 7sc ^ 1, as the estimated zone boundary is in this region. The 
energy of the scattered photon is given hy E — 7sc(l + Psc cosi?) bme(? where •§ is the scattering angle in the electron 
rest frame. The mean expectation for cos'd is zero and the average photon energy after scattering is %cbmeC^. The 
photon direction in the lab frame is nearly aligned with the electron velocity, within an angle i? ~ 7s^^j and hence the 
photon is emitted nearly parallel to the local magnetic field. The photon propagates along the straight line through 
the magnetosphere (the bending of photon trajectory by neutron-star gravity is neglected). We adopt the following 
definition: point P belongs to the radiative zone if a photon (in the || polarization state) emitted parallel to B with 
energy b m^c^ does not convert to a pair and escapes the magnetosphere. 

We have numerically calculated the boundary of the radiative zone for dipole magnetospheres around neutron stars 
with various magnetic moments and temperatures. The boundary is a closed surface around the star, which may 
be described in spherical coordinates by a function ro(0). Instead of r and 9 we find it more convenient to use the 
coordinates i?max and 6, where i?max labels the field lines and the value of b is used as a coordinate running along 
the field line. Then the boundary is described by a function 6o(-Rmax)- The radiative zone only exists on sufficiently 
extended field lines. 



The numerically calculated 60 for i?„iax > ^0 is well approximated by 

60 (-Rmax) = 0.265 ^ A — ^) ^™-l) , (B4) 

^ ' Vl033GcmV lo.SkeV; \ ' ^ ' 



X 0.31 / , \ -0.14 

fl \ kl 



where 

i?* = 5 X 10^ ( ^ ) ( ) cm. (B5) 

VlO^^GcmV V0-5keV; ^ ' 

The corresponding location of the boundary in spherical coordinates, ro(0), is determined by the relations i?max = 
r/ siv?9 and B = (/t/r'^)(l + 3cos^6')^/^. The exact boundary location is reproduced by Equation (IB4p with the 
accuracy of a few percent in a broad range of relevant T (0.2 < kT < 1 keV) and A > 3 x 10^^ G cm . Note that 
i?* appearing in Equation (jB4l) is slightly smaller than Rq. For the models shown in this paper, fi = 5 x 10^^ G cm'^, 
kT = 0.5 keV, and R^ « 0.9i?o- 

We have also calculated the free path of scattered photons to absorption, Zabs, in the pair-loading zone b > bo. It 
depends on the scattering location and the scattering angle. Practically all photons scattered in the zone b > bo have 
lahs/f 1, i-e. they immediately convert to pairs after scattering. At the boundary b = bo, there is a steep 
transition to the escape regime /abs — 00. The small Zabs in the pair-loading zone implies that photons scattered on a 
field line i?max convert to pairs on almost the same field line. For the sample models presented in this paper, we found 
that the typical shift in i?niax before conversion is 5i?max/^max ~ 0.1 for photons scattered at & = &o; ^Rmax quickly 
decreases to essentially zero for photons scattered at & > &o- This fact simplifies calculations in Section 3, as one can 
assume that the processing of the outflow energy into pairs occurs independently on different field lines. 

Note that bo significantly increases near the magnetic axis, where u = R/Rmax 0; sufficiently close to the axis, 
the pair loading zone disappears. This fact is, however, of little importance, as the maximum electric current, which 
corresponds to twist amplitude V' ~ 1: is small near the axis (/ oc u^, see Equation |C7j ). and its contribution to the 
observed luminosity is negligible. The luminosity peaks as at larger u; in this region &o ^ 1/4. 

EMISSIVITY OF THE J-BUNDLE 

Section 3 described how the net emission from the j-bundle is obtained by integrating along thin magnetic tubes 
and then taking the sum over the tubes. It is also useful to have an explicit expression for the spatial distribution 
of emission (e.g. in spherical coordinates) for a concrete magnetospheric configuration. Below we write down the 
general expression for the emissivity per unit volume and then specialize to approximately dipole (moderately twisted) 
magnetospheres. 

First note that a short piece dl of a thin tube with magnetic fiux dF has volume 

dV=^ dl, (CI) 
B 
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where we used the fact that the perpendicular cross section of the flux tube equals B ^dF ■ Then, using dL/dF = 
X^dl/dF (Equation [27]) and dL/dF db — b^^dL/dF (Equation [M]). we find the produced luminosity per unit volume 



dV 2bo 



dl_ dB_ 

dF ~dl 



(C2) 



Here the factor of 1/2 takes into account that db corresponds to two pieces dl of the tube, as the tube has two "legs" 
(two footpoints); $ is the voltage between the footpoints, i.e. the net voltage along the tube. Equation (IC2[) shows that 
dL I dV depends on the distribution of electric current among the magnetic field lines in the j-bundle; this distribution 
is given by 

dF 47ri?,nax' ^ ^ 

where j is the current density. In Equation (jC3[) we used the Maxwell equation V x B = (47r/c)j with neglected 
displacement current, and substituted V x B w tpB / Ry^i^y,. This relation is derived below for the twisted dipole 
magnetosphere; it remains approximately valid for more general configurations. Substitution of Equation (IC3[) to 
Equation (IC2p gives 



dL X* 
dV "1^ 



dB 



dl 



8ttR 



(C4) 



where i?inax is the apex radius of the magnetic field line passing through dV. 

Consider a dipole magnetosphere with a moderate axisymmetric twist. Its magnetic field is well approximated as 
the sum of a dipole poloidal field and a toroidal field B^ (see B09). The poloidal fiux function is that of a normal 
dipole, 

/ ^ - sin^ 6 2nu 



It is also convenient to use the dimensionless fiux function 

_ rji_ R 

27r/i i?nia: 



(C6) 



where R is the radius of the neutron star and i?max is the apex radius of the magnetic field line. Variable u < 1, like 
F, can be used to label the magnetic fiux surfaces or magnetic field lines in the axisymmetric magnetosphere. Suppose 
the region it < is activated, i.e. carries electric currents, and has approximately uniform twist amplitude ^ « const. 
The electric current fiowing between the magnetic axis and a fiux surface u is given by (B09), 

It gives dl/du and dl/dF that confirms Equation (|C3p . One can also derive an explicit expression for dB/dl that 
appears in Equation (|C2p . For the approximately dipole field, B = /tr^^(l + Scos^ 6*)^/^ and 

cLB _ 3/i (3 + 5 cos^ 9) cos 6> 

"dT ^ r4(l+4cos2 6) ■ ^ ' 
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